Gravitational potential evolution in Unified Dark Matter Scalar Field Cosmologies: 

an analytical approach 



O 



Oliver F. Piattella^'^] and Daniele Bertacca^'^-^l 

^ Departamento de Fisica, Universidade Federal do Espirito Santo, 

avenida Ferrari 514, 29075-910 Vitoria, Espirito Santo, Brazil 
^INFN sezione di Milano, Via Celoria 16, 1-20133 Milano, Italy 
^ Dipartimento di Fisica Galileo Galilei, Universita di Padova , via F. Marzolo, 8 F35131 Padova, Italy 
*'INFN sezione di Padova, via F. Marzolo, 8 1-35131 Padova, Italy and 
''Institute of Gosmology & Gravitation, University of Portsmouth, 
Dennis Sciama Building, Portsmouth, POl 3FX, United Kingdom 

We investigate the time evolution of the gravitational potential $ for a special class of non- 
adiabatic Unified Dark matter Models described by scalar field lagrangians. These models predict 
the same background evolution as in the ACDM and possess a non-vanishing speed of sound. We 
provide a very accurate approximation of $, valid after the recombination epoch, in the form of a 
Bessel function of the first kind. This approximation may be useful for a future deeper analysis of 
Unified Dark Matter scalar field models. 
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I. INTRODUCTION 

In the last decade the ACDM model has emerged as the 
concordance model of our universe. It assumes General 
Relativity (GR) as the underlying theory of gravity, and 
two unknown components dominating the late-times dy- 
namics: Cold Dark Matter (CDM), responsible for struc- 
ture formation, and a cosmological constant A making 
up the balance for a spatially flat universe and driving 
the cosmic acceleration Issues (still open) with A, i.e. 
the cosmological constant problem, prompted the cosmol- 
ogists' community to consider alternatives, mainly in the 
form of a dynamic component dubbed Dark Energy (DE) 
(see, for example, However, it should be recognized 
that, while some form of CDM is independently expected 
to exist within modifications or extensions of the Stan- 
dard Model of high energy physics (for example, see 
the really compelling reason to postulate DE has been 
the acceleration in the cosmic expansion [l| . It is mainly 
for this reason that it is worth investigating the hypothe- 
sis of CDM and DE being two aspects of a single Unified 
Dark Matter (UDM) component. 

A large variety of UDM models, mainly based on adi- 
abatic fluids or on scalar field Lagrangians have been 
investigated in the literature, see jj] for an up-to-date 
review. The most important feature appears to be the 
existence of pressure perturbations in the rest frame of 
the UDM fluid. The latter correspond to a Jeans length 
below which the growth of small density inhomogeneities 
is impeded and the evolution of the gravitational poten- 
tial is characterized by an oscillatory and decaying be- 
havior. Therefore, the general lesson that may be drawn 
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is that a viable UDM model should be characterized by 
a small speed of sound (typically < 10~'*c i6|]), so that 
structure formation would not be hampered and the ISW 
effect signal would be compatible with CMB observation 
(see also [1^). 

Among those approaches which describe UDM models 
as scalar fields, we focus on the one presented in fl^ . 
Here, the authors devise a reconstruction technique for 
models whose speed of sound is so small that the scalar 
field can cluster and, at the same time, the background 
expansion is identical to the one of the ACDM. This class 
of models has recently been investigated in by means 
of the cross-correlation between the ISW effect signal and 
the large scale structure (LSS) distribution and in 6j by 
means of parameters forecasts for future 3D cosmic shear 
surveys (see also [H]). In this brief report, we investigate 
in detail the evolution of the gravitational potential for 
this class of models. Through some simple mathemat- 
ical manipulations we show how to cast the usual time 
evolution equation in a form which can be very reason- 
ably approximated by a Bessel type equation, after the 
recombination epoch. The resulting solution provides a 
fitting function in excellent agreement with the numerical 
results. 

Throughout the paper we use SttG = c? = 1 units and 
the (+,—,—,—) signature for the metric. 



II. THE MODEL AND THE BASIC EQUATIONS 

Assuming a fiat Friedmann-Lemaitre-Robertson- 
Walker (FLRW) background metric with scale factor 
aij]) and 77 the conformal time, the authors of [l2| 
introduced a new class of non-adiabatic UDM scalar 
field models which, by allowing a pressure equal to —A 
on cosmological scales, reproduce the same background 
expansion as the ACDM one. When the energy density 
of radiation becomes negligible, and disregarding also 
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the baryonic component, the background evolution of 
the universe is completely described by 



(1) 



where H = a' /a is the conformal time Hubble parame- 
ter and the prime denotes derivation with respect to the 
conformal time; Hq is the Hubble constant and f^Ao and 
^mo = 1 — ^^Ao are to be interpreted as the "cosmological 
constant" and "dark matter" density parameters, respec- 
tively. From the WMAP7 data: SIaq = 0.734± 0.029 'H. 
Moreover, the model is characterized by a speed of sound 
[TgI [17| which has the following form : 



III. THE EQUATION FOR T 

Inserting Eq. ^ into Eq. (|4]) and using the equation 
for "I>A one obtains 



T' + clk^T = 



(7) 



Being such equation not analytically solvable, in [13| the 
authors employed the following approximation for T: 



cs{v')dv' ■ (8) 
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where Cj^ 



(a 



limit and v = 



IS its asymptotic 
i 0.362. 

Consider small inhomogeneities of the scalar field and 
write the perturbed FLRW metric in the longitudinal 
gauge jlT]: 

ds^ = a{r]f [(1 + 2$) dr]^ - (1 - 2$) S.jdx'dx^ , (3) 

where $ = $(x, ry) is the gravitational potential. Con- 
sidering plane- wave perturbations $ cx exp (ik • x) of co- 
moving wave-number k = |k| the evolution equation of 
$ becomes 



3-H$' 



^fc2$ + (2-H' + -H^)'I> = 0, (4) 



where Cg is given in Eq. Note that the gravitational 
potential evolves in the same way as in a ACDM universe 
for those modes with wavelengths larger than the sound 
horizon and/or when Cg — >■ 0. In the same spirit of 
we express the gravitational potential $ {k; rj) of Eq. 
as follows: 



$ (fc; ry) = $A (fc; r/) T (fc; ?7) 



(5) 



namely as a transfer function T (k ; r/) applied to $a, 
where $a is the ACDM gravitational potential, which 
solves Eq. Q with Cg = 0. In particular, we impose that 
T{k; r]) = 1 and dT{k; 'r])/dr] = for 77 < ?7rec, where ryroc 
is some epoch when the universe is matter dominated and 
the radiation component is negligible (usually later than 
the recombination epoch). Before this time, the sound 
speed must be very close to zero to allow for structure 
formation. 

The knowledge of T (fc ; rj) is useful because it allows to 
directly translate informations, such as the total matter 
power spectrum, from the ACDM to the UDM models 
of [T^]. In order to compare the predictions of the mod- 
els with the observational data, one should express the 
density contrast as a function of the gravitational poten- 
tial. For example, for scales smaller than the cosmologi- 
cal horizon and z < Zj-cc one has 



Sum [kwi^)] = 



-2fc2$A [k:'q{z)]T[k:Tj{z)] 

3i/2f^mO (1 + Z) 



(6) 



where (5dm is the UDM density contrast [l^, [13 ■ 



The integral defining Aiji) is not solvable, but since we 
are going to consider < 10"^ 0j[i11j it is convenient 
to expand Cg in Taylor series near Coo = 0: 



l + iy{l + zy 



(9) 



Retaining just the first order of this expansion, plugging 
it into the formula for A{ri) in Eq. ^ and changing the 
integration variable to the scale factor a, one finds 



A: 



da 



7-3 



da- 



(10) 

For simplicity, choosing as lower integration limit a = 0, 
the integration can be performed analytically and the 
result is: 
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As graphically showed in [ijj, Eq. ([8]) together with 
Eq. (jlip gives quite an accurate approximation of T for 
this class of UDM models. In the following, we give a 
motivation for such good agreement and our discussion 
shall allow to find an even more accurate fitting function. 
In the next section, we analyze the term $^/$a, which 
plays an important role in Eq. ©. 



IV. THE ACDM GRAVITATIONAL POTENTIAL 

In terms of the scale factor a, the evolution equation 
for $A becomes 



cP$A 



n' 



aU^ 



d$A 

da 



1 



2%' 
am^ 



*A = . (12) 



Considering a -> ^ = {1- ^/oF/l7+l) /2, we cast Eq. ([12 
in the form of a hypergeometric equation: 



e(C-l)' 
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d^ 



(2^-1 



(13) 
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Choosing $A(fc;arec) = 1 and d<^f^/da\^ = as nor- 
malized initial conditions, and applying a suitable Kum- 
mer transformation [l^ one can cast the solution in the 
following form: 

$A(fc;a) = 4>(fc;0)T^(fc)2Fi , (14) 

where $ (fc; 0) is the primordial gravitational potential at 
large scales, set during inflation [l^, and T^{k) is the 
matter transfer function up to recombination suggested, 
for example, by or \21] (for a more general and deeper 
analysis of the ACDM model, see also [22j). From ([T4|. 
we are able to compute the term ^^/^a hi Eq. ([7]). 



In Figs. [2] and |3l the fitting function (ITTl) is compared 
with the fitting function of Eq. , the same adopted in 
[l3l |. For both the cases, the functional form of yl(r/) is 
the approximated one given in Eq. dTTT) . We choose, as 
a representative example, the case — 10^^, for which 
the oscillatory behavior is most pronounced. 

The reason for the success of 5 = 9/4 is the following: 
on average, in the interval a € (10"'^, 1) we do have b sa 
9/4. Indeed, computing the average of b with respect to 
a or g, we approximatively obtain {b)a ~ {b)g « 2.30, 
where 




V. THE FITTING FUNCTION 

Introducing a convenient new variable g{ri, k) = kA{ri), 
Eq. (O becomes 

g^^+gbig,7j)^+g^T^0, (15) 



where 



u \ - 9 9 ^ g /„*A 
9 9 \ 



(16) 
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FIG. 1: Plot of b as function of the scale factor a and for 
= 0.362. 



In Fig. [T]we plot b as function of the scale factor and for 
1/ = 0.362. Notice that 6 w 9/4 for most of the evolution, 
i.e. up to a w 0.5. Then, it varies from 9/4 (2.25) to 
2.47 at a = 1. By approximating b — 9/4 we make a 
relative error of about 9% in a = 1. Remarkably, the 
choice 5 = 9/4 renders Eq. p3|) a Bessel-type equation 
llSl. whose solution is 



T(5) = 2'V8r(H)5-tj,(g) 



(17) 



FIG. 2: (Color online). Upper panel: comparisons for = 
10~^ and k = 0.1 h Mpc~^ of the numerical result (black 
solid line) for T, the fitting function ((8]) (red dashed line) 
and the fitting function H17p (blue dash-dotted line). Lower 
panel: absolute differences between the numerical result for T 
and the fitting function ((8| (red dashed line) and the fitting 
function (|17|) (blue dash-dotted line). 



In conclusion, from Eq. p7)) and the definition of 
Sbm {k;z), in Fig. |3]we show the power spectrum P{k) 
of the UDM energy density component that clusters: 

{Sum (k; z = 0) ,5dm (k'; ^ - 0)) = (2^)35(k + k')P{k) , 



(19) 
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FIG. 3: (Color online). Upper panel: comparisons for — 
W^^ and a = 1 of the numerical result (black solid line) for 
T, the fitting function ^ (red dashed line) and the fitting 
function p7l) (blue dash-dotted line). Lower panel: absolute 
differences between the numerical result for T and the fitting 
function ((8| (red dashed line) and the fitting function p7p 
(blue dash-dotted line). 



see also Il3ll. 



VI. CONCLUSIONS 



obtaining a new fitting function, see Eq. p7)) . that de- 
scribes with great accuracy the evolution of T (fc ; 77) after 
tlie recombination epoch. This fitting function provides 
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FIG. 4: (Color onhne). Power spectrum P{k) of the UDM 
clustering energy density component 5dm. From left to right, 
= 10"^ 10"*, 10"'^, 10"'' (green solid, blue dashed, red 
dot-dashed, black dotted lines, respectively). 



a useful tool for testing this class of models, especially in 
the light of new data on the cosmic background radiation 
and the weak gravitational lensing coming from, among 
the others, the Planck collaboration [23] and the Euclid 
project [2I I. 

A future development of our work would be to inves- 
tigate in more detail the evolution of the gravitational 
potential in UDM models characterized by a fast transi- 
tion [13, Eli- In the adiabatic case, the speed of sound is 
a function highly peaked in the moment of the transition, 
being exponentially vanishing elsewhere. We expect the 
analytic approximate treatment to be more complicated 
than the one here presented, because it would involve 
mathematical techniques for dealing with rapidly varying 
coefficient functions in a second order differential equa- 
tion. 



We have investigated in detail the evolution of the 
gravitational potential for a class of non-adiabatic Uni- 
fied Dark Matter models which provide a background 
evolution similar to the ACDM one and possess a non- 
vanishing speed of sound. Such class of models has been 
investigated in [l^ where a convenient fitting function 
for the gravitational potential, in the form of a spheri- 
cal Bessel function of order zero, has been employed. In 
particular, we have analyzed the evolution equation of a 
suitable transfer function T (fc ;77), defined in ([5]), intro- 
ducing the new variable g — kA{ri), cf. Eq. ([15]), and 
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